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thermal energy storage with phase change materials, the volume of literature begins to be significant. This 
investigation carries out a compilation of information on two latent thermal fluids: PCM emulsions and 
microencapsulated PCM slurries (mPCM slurries). This review presents tables containing information on 
the different PCM emulsions and mPCM slurries studied by different researchers, as well as commercially 
available products. Thermophysical and rheological properties are analyzed, making a special effort to 
analyze heat transfer phenomena, concluding with the enumeration of application examples available in 
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1. Introduction 


There is a clear necessity to rationalize the use of energy in 
modern society. One of the main contributions to the rational use 
of energy is thermal energy storage. There is a great technical 
potential in the use of thermal energy storage substituting the com- 
bustion of fossil fuels, allowing for a greater use of energies which 
do not present constant supply, such as residual heat or renew- 
able energies. Energy sources can be more efficient using short- 
and long-term storage. Thermal energy storage can also be used to 
shave or eliminate the peak of electricity demand, including more 
expensive electricity that is generated during on-peak hours. 

Thermal energy storage based on the use of the solid-liquid 
phase change of materials (PCM) is a technique that has been 
recently raising a greater practical interest. Such an interest is moti- 
vated by the great thermal energy storage density per unit volume 
of PCM materials in a reduced temperature range, and also by the 
constant incorporation of new materials with very different proper- 
ties and phase change temperature intervals. The interest in phase 
change materials is evident when considering thermal energy stor- 
age systems with PCM, which have been developed for different 
applications [1]. 

Recently, a new technique has been proposed to use phase 
change materials in thermal storage systems, heat exchangers and 
thermal control systems. This new technique consists of forming a 
two-phase fluid, from the mixture of a fluid, such as water, and 
a phase change material, such as paraffin, resulting in a latent 
heat storage fluid. Inaba [2] classifies thermal fluids, describing the 
main characteristics and applications. Among the latent thermal 


fluids, five types of fluids are mentioned: 1) ice slurries; 2) phase 
change material microemulsions, in which the PCM is dispersed 
in water through an emulsifying agent; 3) microencapsulated PCM 
slurries, where the PCM is microencapsulated in a polymeric cap- 
sule and dispersed in water; 4) clathrate hydrate PCM slurries, 
where the clathrate hydrates are composed of water molecules 
(host molecule) forming a weaved structure where the molecules 
of the other substance (guest molecule) are accommodated, consti- 
tuting a special molecular structure where the heat associated with 
the chemical reaction of formation and dissociation of clathrate 
hydrate is greater than that of ice melting and 5) shape-stabilized 
PCM slurries (ssPCM slurries), based on ssPCM, these can consist 
of paraffin infiltrated in high density polyethylene, with a melting 
temperature higher than of the paraffin. In this way the paraffin is 
retained inside the structure of high density polyethylene, avoiding 
the leak of the PCM. Fig. 1 shows a schematic draw of the different 
types of PCM slurries [3], although the case of clathrate hydrate PCM 
slurries does not appear drawn. In 2010, Zhang et al. [4] published 
a review, focused on two latent thermal fluid types: mPCM slurries 
and clathrate hydrate slurries, putting special emphasis on ther- 
mal properties and applications. The review presented herein aims 
at completing the work of Zhang et al. [4] regarding mPCM slur- 
ries and going a step further by analyzing the information available 
in literature on PCM emulsions and by studying the heat transfer 
phenomenon. 

As main issues to be tackled, some studies inform that in the 
case of mPCM slurries it is particularly difficult to maintain a stable 
homogeneous flow ifthe particles are not processed with very small 
size and high flexibility. Besides, the PCM capsule entails an extra 
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Fig. 1. Schematic drawing of the different types of PCM slurries [3]. 
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cost. The capsule prevents the PCM in continuous phase from leak- 
ing, which in that case could solidify in ducts and cause clogging. 
It is important that the capsules are sufficiently resistant against 
the stress produced by the pumps. In the case of emulsions, pre- 
vious experiments indicated that it is difficult to maintain a stable 
emulsion above melting temperature, as instabilities could appear 
during phase change [5]. Stratification problems will appear, as 
the paraffin droplets will form greater droplets and finally a PCM 
layer will float in the upper part of the storage system, due to the 
difference of densities [3]. 

These new fluids offer many advantages and can be used either 
as thermal storage materials or heat transfer fluids [6], due to: 1) 
high storage capacity during phase change, 2) possibility to use the 
same medium either to transport or store energy, as these slur- 
ries are pumpable (reducing in this way heat transfer losses), 3) 
heat transfer at an approximately constant temperature, 4) high 
heat transfer rate due to the elevated ratio surface/volume, 5) 
lower pumping power, as a consequence of the reduction in mass 
flow due to higher heat capacity, 6) a better cooling performance 
than conventional heat transfer fluids, due to the decrease in fluid 
temperature as a consequence of higher heat capacity. As addi- 
tional advantages, these novel fluids present a thermal energy 
storage density that results advantageous in comparison to con- 
ventional systems of sensible heat storage in water and can be 
competitive against macroencapsulated PCM tanks. According to 
the consulted literature, macroencapsulated PCM tanks present 
porosity values between approximately 0.37 and 0.47 [7-9]. This 
means that the PCM volume together with the container system 
represents approximately 53-63% of the volume of the thermal 
energy storage system. Therefore slurries with PCM concentra- 
tions volumes of approximately 30-40% could result competitive 
as thermal energy storage systems. Besides, response time will 
be shorter using these PCM emulsions or mPCM slurries as stor- 
age material than with macroencapsulated PCM. The tanks will 
be simpler as there is no need to macroencapsulate, and conven- 
tional tanks can be used. Nevertheless, currently there is a lack of 
experience on its technical viability. In order to be advantageous, 
these latent fluids must meet the following requirements [10]: 1) 
high heat capacity, 2) phase change temperature range matching 
the application, 3) low subcooling, 4) high heat transfer rate, 5) 
pumpable, low pressure drop in pump systems, 6) stable over along 
term storage, and 7) stable to thermal-mechanical loads in pump 
systems. 

This article encompasses a review of scientific literature as well 
as a technical and commercial compilation of PCM emulsions and 
mPCM slurries regarding their use as heat transfer fluids and ther- 
mal storage materials. This article completes the review by Zhang 
et al. [4] with the collection of additional information on mPCM 
slurries and by considering PCM emulsions. Several tables present 
information on PCM emulsions and mPCM slurries (PCM micro- 
capsules) studied in scientific literature, as well as commercially 
available products, presenting the main characteristics of the mate- 
rials. This review goes a step further by making an exhaustive 
analysis on heat transfer phenomena in these fluids, concluding 
with the factors affecting the phenomena. The review is structured 
in five greater sections: a first section that considers the elaboration 
of emulsions and the different fabrication techniques for microcap- 
sules; a second section that presents tables containing compilation 
of PCM microcapsules and emulsions, both in research stage and 
commercially available; a third section that analyzes the thermo- 
physical and rheological properties of PCM emulsions and mPCM 
slurries, also considering the subcooling problem and stability of 
fluids; a fourth section that deepens knowledge on heat trans- 
fer phenomena; and finally, a fifth section that enumerates and 
describes application examples where these new fluids have been 
employed. 


2. Manufacture process 
2.1. Elaboration of emulsions 


As previously explained, emulsions are systems formed by two 
unmiscible liquid phases, one of which is dispersed in the other 
in the form of droplets through substances that modify interfacial 
activity, denominated surfactants. In the case of emulsions, spe- 
cial attention must be given to viscosity and stability. The viscosity 
of an emulsion will increase when thickening agents are aggre- 
gated, avoiding unstabilization processes by increasing the disperse 
phase proportion and reducing the size of the emulsion particles. 
Emulsifying methods can be classified in two groups: high and low 
energy methods [11]. High energy methods produce coarse emul- 
sions (5-100 pm) through a high degree of shearing generated by 
helix or turbine mixers. A greater reduction in the droplet size 
requires a great amount of mechanical energy. The formation of 
nanoemulsions using high energy methods is usually very expen- 
sive, as it requires the use of high pressure homogenizers, ultrasonic 
generators, microfluidizers, etc. The low energy methods are those 
where the energy required to form the emulsion originates from 
transitions or phase changes occurring during the production of 
the emulsions. 


2.2. Fabrication of PCM microcapsules 


The microencapsulation technique has been widely used in the 
pharmaceutical and chemical engineering fields. In recent years 
this technology has reached the field of phase change materials 
in order to improve their behavior. Microcapsules are small par- 
ticles of material coated by another material, the second material 
forming a thin film over the first one, isolating and protecting from 
the environment. Size range of microcapsules is quite wide, with 
diameter slurry between 2 and 2000 pm. The size of the shells oscil- 
lates between 0.5 and 150 pm thickness and the core constitutes 
between 20 and 95% of the total mass. There are many different 
fabrication techniques for microcapsules, and the choice is made 
depending on the characteristics of the active material to be encap- 
sulated and the type of polymeric material used as shell [12]. 

Ina first stage, a dispersion is formed with the active material. In 
a second stage, this dispersion is transformed and finally, stabiliza- 
tion and solidification techniques are applied to the shells before 
separation. Emulsifying techniques are, in principle, the most sim- 
ple and extended for the obtaining of microcapsules. 

There are several methods that can be used to produce micro- 
capsules. Depending on the nature of the process, there are 
physical, physical-chemical, and chemical processes. The most 
utilized techniques in literature for PCM microencapsulation 
according to the aforementioned classification are: for the first 
group, the spray drying technique; for the second group, coacerva- 
tion; and for the third group, in situ and interfacial polymerization. 
These techniques are explained in the following paragraph. 

The spray-drying technique is based on the preparation of an 
emulsion, dispersing the material of the core (in this case the PCM) 
in a concentrated solution of the material forming the capsule, 
until the desired particle size is obtained. This emulsion is pulver- 
ized into droplets, which will be next dried when makes contact 
with a hot medium, evaporating the present water instantly, and 
consequently allowing the active material to be trapped inside a 
film of encapsulating material. In the coacervation technique, the 
polymeric solute is separated in the form of small liquid droplets 
(forming the coacervate) and deposited around the insoluble par- 
ticles dispersed into a liquid. These droplets slowly unite and form 
a continuous cover around the core. In order to obtain a longer 
lifetime of microcapsules, coacervation can be accomplished in 
two steps, adding the polymer twice. Through this mechanism, a 
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Fig. 2. Morphology of microcapsules through a SEM microscope [13]. 


smaller thickness of the microencapsulate shell can be obtained 
with the same weight index, during a longer time so that the poly- 
mer regulates its molecules with the PCM core. Compactability and 
impermeability are improved (a lower speed of polymer deposi- 
tion increases impermeability), proportioning greater stability to 
the microcapsules, conserving size and spherical form. The texture 
is smoother and the spherical form is more regular, compared to 
one-step coacervation, where microcapsules with many protru- 
sions, rougher, coarser and more porous are obtained [13]. The 
morphology can be observed in Fig. 2. In the interfacial polymeriza- 
tion technique, the polymer constituting the microcapsule cover is 
formed from two monomers in separate phases (an aqueous phase 
and an organic phase) that react in the interface of both phases, 
at the moment of microencapsulation. In the in situ polymeriza- 
tion, firstly a PCM emulsion is prepared and then the synthesis of 
the prepolymer solution is carried out through the mixture of two 
polymers, which will form the cover, and water. This prepolymer is 
added to the emulsion in the form of droplets, while the emulsion is 
agitated during a specific time. The emulsion is cooled and filtered, 
obtaining the microcapsules, which have to be dried. 

According to the study of Yang et al. [14], the polymers that 
present the best characteristics at the time of microencapsula- 
tion are the flexible plastics, such as the Polymethyl methacrylate 
(PMMA) and polyethyl methacrylate (PEMA). With these plastics, 
a higher value of melting and solidification latent heat is achieved. 
Besides, the viscosity of the slurry is not significantly affected by 
the material of the microcapsule shell. The main manufacturer of 


PCM microcapsules and mPCM slurries is BASF. BASF manufactures 
its microcapsules from a polymerization process of a PCM emul- 
sion. This microencapsulation process can be observed in Fig. 3 [15]. 
Firstly a paraffinic emulsion is elaborated in hot water with agita- 
tion and detergents. Then the monomers are added and the growing 
process of the polymeric capsule begins in the interface of paraf- 
fin and water. Finally a highly-reticulated polymer forms a dense 
capsule over each paraffin droplet. This process results in a liquid 
product, in which the capsules are dispersed in water. The micro- 
capsules are obtained after the dispersion passes through a drying 
tower, as shown in Fig. 4. The diameter of the dispersion microcap- 
sules is between the range 2-20 um. Nevertheless, after the drying 
process, secondary particles with a larger size are created. 


3. Compilation table of PCM emulsions and mPCM slurries 


Table 1 shows information on the PCM microcapsules and 
mPCM slurries studied in literature, together with some thermo- 
physical properties, as well as additional interesting information 
given by the authors. Table 2 shows a compilation of PCM emulsions 
studied in literature. 

At the moment there are few PCM manufacturers approaching 
these techniques. Table 3 collects data on commercially available 
PCM emulsions, mPCM slurries and PCM microcapsules up to date. 

As is known, inorganic PCM, usually salt hydrates, present 
certain advantages when compared to organic PCM, such as 
a higher energy storage density, higher thermal conductivity, 


Table 1 


PCM microcapsules and mPCM slurries studied in literature. 


Ref. Microencapsulation Core material Shell material Nucleation agent Size distribution Melting temperature Phase change enthalpy 
process 
[13] Two step coacervation n-Octadecane Melamine formaldehyde 2 um 
[14] In situ polymerization Tetradecane 1) PVAc (polyvinyl acetate) 2) 5-40 um 1) No manisfested 2)~0 
2) PS (poliestirene) 3) 5-30 um phase change 3) 66.26 kJ/kg 
3) Polymethy! methacrylate 4) 5-30 pm 2) 2.06°C 4) 80.62 kJ/kg (melting 
(PMMA) 3)5.97°C enthalpy) 
4) Polyethyl methacrylate (PEMA) 4) 5.68°C 
Tetradecane content of capsule: 
40% 
[16 Lauryl alcohol Melamine formaldehyde 5-10 um 24°C 
[17 n-Octadecane 1-2 pm 
Ø average: 1.3 ym 
[18 Interfacial polymerization n-Octadecane Poliurethane 5-10 um 
[19 Interfacial polymerization n-Octadecane Urea Melamine formaldehyde 1-Tetradecanol 160 kJ/kg 
n-Nonadecane paraffin 
n-Eicosane 1-Octadecanol 
[20 In situ polymerization n-Octadecane Urea Melamine formaldehyde 0.2-5.6 um 
[21 In situ polymerization Paraffin Urea-formaldehyde 5-20 um (depending ~54°C 157.5 kJ/kg (76.9% core 
on the amount of content) (melting enthalpy) 
emulsifier) 
[22 Polymerization of emulsion Docosane PMMA (polymethyl methacrylate) 0.14-0.466 wm 41°C 54.6 kJ/kg 
Ø average: 0.16 um 
[23 In situ polymerization of a n-Octadecane Poliestirene 50-200 nm 124.4kJ/kg 
nanoemulsion assisted by Ø average: 124nm 
ultrasound 
[24 In situ polymerization n-Octadecane Melamine formaldehyde 30.5°C 170 kJ/kg (melting enthalpy) 
[25 Polymerization n-Octacosane PMMA (polymethyl methacrylate) 0.15-0.33 um Ø 50.6°C 86.4kJ/kg 
average: 0.25 um 
[26 Coacervation Paraffin wax (Merck) 145-240 kJ/kg (depends on the 
Spray-drying ratio core/shell and on the 
method) 
[27 In situ polymerization n-Octadecane Melamine formaldehyde Ø average: 2.2 uym 40.6°C 144 kJ/kg 
[28 Coacervation 1) 99.8% Gelatin Silicon fume 1) 90-150 pm 3) 202.1 kJ/kg (melting 
Tetradecane + 0.2% Tetradecanol Ø average: 100 um enthalpy) 
silica fume 2) 70-260 um 
2) 98% tetradecane + 2% Ø average: 145 ym 
tetradecanol 3) 2-10 um 
3) 94% tetradecane + 6% Ø average: 4.4 ym 
tetradecanol 
[29 n-Octadecane 1-5 um Melting range: 147.1 kJ/kg (melting enthalpy) 
Average diameter 24-29°C 
4.97 ym 
[30 1-Bromohexadecane Amino plastic Average diameter: 14.3°C (melting 5 wt% 6.5 kJ/kg 
Ratio core-shell: 7 8.2 um starting temperature) 10 wt% 13 kJ/kg 
Capsule thickness: 0.3 wm 15.8 wt% 20.5 kJ/kg (melting 
enthalpy) 
[31 RT6 Melting range: 4-6.8 °C For a 45 wt% concentration: 
55 kJ/kg (melting enthalpy) 
[32 DPNT06-0182 (Ciba 10-100 um ~35°C 96.968 kJ/kg 
Specialty Chemicals) 10% 8.074 kJ/kg 
25% 15.194kJ/kg 
35% 37.213 kJ/kg (melting 
enthalpy) 
Micronal DS 5008X 1-20 um ~29°C 102.008 kJ/kg 


(BASF) 
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Table 2 


PCM emulsions studied in literature. 


Ref. PCM Nucleation agent Surfactant Emulsifying method Size distribution Melting temperature Phase change enthalpy 
[5] Mixture of n-alkanes Non-ionic surfactant Ultrasonic generator 2 um 9.5°C 50% 78.9 kJ/kg (melting 
Power: 500W enthalpy) 
Frequency: 20 kHz 
[11 20% Tetradecane 6% surfactant (67.7% Phase inversion 200-250 nm 43 kJ/kg (melting 
Tween60, 32.3% temperature method enthalpy) 
Span60) 
[14 Tetradecane Average diameter 10% 5.06°C 10% 18.5 kJ/kg 
10% 18.18 wm 20% 5.84°C 20% 112.3 kJ/kg 
20% 18.56 um 30% 5.84°C 30% 150.8 kJ/kg 
30% 16.05 um (melting enthalpy) 
[33 Hexadecane <0.1mm 16.5°C 
[34 C22H46 10-40 pm 
[35 Mixture of hexadecane 
and tetradecane 70/30 
[36 30% RT6, RT10, RT20 2.5% paraffin with a 1.5% alcohol ethoxylate For a temperature 
melting temperature of range of 6°C and 
50°C 30 wt% concentration: 
RT6=75 kJ/kg 
RT10=50kJ/kg 
RT20 = 44 kJ/kg (total 
capacity of storage) 
[37] Hexadecane - SDS - Disperser 0.05-30 pm 
- Tween - Ultrasonic generator (depending on the 
method) 
[38] - tetradecane - SDS - Ultrasonic generator 0.05-30 pm (for the 15.4-17.4°C (for the 
- hexadecane - Tween40 - Rotor-stator system hexadecane emulsion hexadecane emulsion) 
- RT20 - Surfactant mixture depending on the 
method) 
[39] Paraffin Phase Incursion Average diameter: 
Method 0.304 um 
[40] 30% RT 10 Melting range: 55 kJ/kg (5-11 °C) (total 


4-11.5°C 


capacity of storage) 
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Fig. 3. Microencapsulation process of BASF [15]. 
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non-flammability, and are cheaper. Nevertheless, inorganic PCM 
also present a series of inconveniences, such as the corrosion prob- 
lem, phase segregation and subcooling. As a consequence of this 
last problem (which in microscopic geometries would be worse 
due to probabilistic reasons), there are still few studies on inor- 
ganic PCM microencapsulation. After systematic reviews, the only 


Table 3 
Commercially available PCM microcapsules. 


manufacturer of inorganic microencapsulated PCM is Capzo, sell- 
ing microcapsules of salt hydrates and the only consulted work is 
that of Salaün et al. [41]. This last work investigates the effects of 
the preparation conditions (solvent evaporation-precipitation pro- 
cess) on the morphology of the capsule and on the efficiency of the 
sodium phosphate dodecahydrate encapsulation. 


Manufacturer Product Type of product PCM Concentration Particle/droplet Melting Latent heat 
size temperature 
BASF DS 5000 mPCM slurry Paraffin 42% 26°C 45 kJ/kg 
DS 5007 mPCM slurry Paraffin 42% 23°C 41 kJ/kg 
DS 5030 mPCM slurry Paraffin 42% 21°C 37 kJ/kg 
DS 5001 Powder Paraffin 26°C 110kJ/kg 
DS 5008 Powder Paraffin 23°C 100 kJ/kg 
DS 5030 Powder Paraffin 21°C 90 kJ/kg 
Microtek Laboratories MPCM-30D Powder n-Decane 17-20 pm —30°C 140-150 kJ/kg 
MPCM-10D Powder n-Dodecane 17-20 pm —9.5°C 150-160 kJ/kg 
MPCM 6D Powder n-Tetradecane 17-20 um 6°C 157-167 kJ/kg 
MPCM 18D Powder n-Hexadecane 17-20 um 18°C 163-173 kJ/kg 
MPCM 28D Powder n-Octadecane 17-20 pm 28°C 180-195 kJ/kg 
MPCM 37D Powder n-Eicosane 17-20 pm 37°C 190-200 kJ/kg 
MPCM 43D Powder Paraffin mixture 17-20 pm 43°C 100-110kJ/kg 
MPCM 52D Powder Paraffin mixture 17-20 um 52°C 120-130 kJ/kg 
Capzo Thermusol HD35SE Powder Salt hydrate 30-40°C 200 kJ/kg 
Thermusol HD60SE Powder Salt hydrate 50-60°C 160 kJ/kg 
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4. Main characteristics of PCM emulsions and mPCM 
slurries 


4.1. Solidification, hysteresis and subcooling 


When different results between the measurements of the cool- 
ing and heating processes are obtained, this phenomenon is called 
hysteresis; nevertheless, hysteresis is commonly understood as a 
material property. The effects (due to the measurement conditions) 
are called apparent hysteresis. In the case of inorganic materials the 
hysteresis phenomenon is much more noticeable than in the case 
of organic materials. 

There are several effects that originate hysteresis due to the 
material, and the most common is subcooling. When a liquid has 
to be cooled to a temperature below its melting temperature so 
that the crystallization process to begin, this process is called sub- 
cooling. Subcooling is always a serious issue in PCM investigation 
and application fields, as the operation temperature range of the 
storage systems will be wider, worsening the energy efficiency of 
the systems. Despite the fact that many PCM do not present the 
subcooling phenomenon in macroscopic geometries, the problem 
occurs in microscopic geometries. The hysteresis and subcooling 
problem can be observed in Fig. 5. 

Phase transformation, in this case solidification, can be 
described as a nucleation process and posterior growth of the stable 
nuclei formed. Depending on the nature of the material and on the 
transformation, two types of nucleation can occur: homogeneous 
and heterogeneous nucleation. Homogeneous nucleation occurs if 
all volumes of the initial phase are identical from energy, chemical, 
and structural viewpoints. This is not the most common process 
as materials inevitably contain defects or non-homogeneities that 
cause local differences in some properties. In this case nucleation 
will occur preferably in those locations where non-homogeneities 
occurred. The process ceases to be random as there are preferen- 
tial nucleation locations and then heterogeneous nucleation takes 
place. In heterogeneous nucleation the formation of nuclei begins 
in preferential locations such as walls or imperfections. The addi- 
tion of impurities such as nucleation elements is a technique used 
to obtain greater homogeneity in the structure of the solidified 
material. The aggregation of particles that act as nucleation agents 
can result in solidification at lower subcooling, avoiding the forma- 
tion and growth of nuclei mainly from the walls of the mold, and 
leading to structures with anisotropic behavior. In his experimen- 
tal work, Huang et al. [36] utilized paraffin as a nucleation agent, 
with a higher melting temperature than the material in solidifi- 
cation, obtaining in this way solid surfaces on which to begin the 
solidification process. Nevertheless, due to the microscopic size of 
particles in the case of PCM emulsions and mPCM slurries, the prob- 


ability that nucleation agents are present to begin the solidification 
process of the PCM is lower, and therefore a lower temperature is 
required for solidification to start. This effect is explained by Giin- 
ther et al. [37] in her article. For a given concentration of seeds, the 
number of active volumes (droplets) by these seeds is dependent on 
the size of the volumes. The remaining inactive volumes (droplets 
without seeds) can only solidify after homogeneous nucleation. The 
fraction of active volumes with respect to the inactive volumes in 
function of the size of the element corresponds to a Poisson distri- 
bution. 

There are many authors that have studied the subcooling prob- 
lem in PCM emulsions and mPCM slurries as well as the way 
to eliminate such problem. The authors of the different studies 
consider the value of subcooling as the difference between the 
peak temperatures in the heating curve and in the cooling curve, 
obtained with the DSC. 

In the case of the paraffin emulsion developed by Huang et al. 
[10], in order to prevent subcooling, a nucleation agent was added 
to the paraffin before the emulsifying process. The nucleation agent 
offered the necessary surface to start nucleation in the interior of 
the paraffin droplets, acting therefore as a nucleation catalyser. The 
emulsion without nucleation agents presented a 7°C subcooling, 
while with the addition of the nucleation agent it was reduced to 
0.1 °C. To study the effect of the droplet size on the nucleation and 
subcooling, Huang et al. [38] prepared different paraffin-in-water 
emulsions with various droplet sizes. They found that besides the 
droplet size, the surfactant of the emulsion had a significant influ- 
ence on the solidification. They also observed that all emulsions 
containing nucleation agents showed a freezing peak spectrum. In 
the second part of this paper [42], Günther et al. discussed possi- 
ble interpretations of the observations according to the nucleation 
theory. They relate the increased subcooling in emulsions to the 
droplet size. The seed deactivation for a given seed concentration 
increases when the droplet size decreases. Regarding to the surfac- 
tants, they point out that their chemical structure could promote 
the nucleation. They also mention that while in larger droplets (low 
curvature), the alignment of the surfactants tails is almost paral- 
lel, offering in this way a pattern for crystallization of the alkane 
crystal, in smaller droplets (high curvature) this alignment is more 
complicated. The freezing peak spectrum is explained with differ- 
ent seed types, as it can be observed in Fig. 6. In a bulk sample, 
nucleation is dominated by primary seeds (dots), at their nucle- 
ation temperature T;. When the volume is divided into droplets, 
secondary seeds (squares) become active at T2. For smaller droplets, 
nucleation on the droplets wall or homogeneus nucleation comes 
into the play at T3. Alvarado et al. [43] studied the subcooling pro- 
cess in n-tetradecane microcapsules, as well as its elimination due 
to the addition of nucleation agents. N-tetradecane microcapsules 
were prepared with gelatine capsules, producing microcapsules 
with an average diameter of 100 um, utilizing two different sub- 
stances and concentrations as nucleation agents (0.2% silica, 2 and 
4% tetradecanol). Silica fume did not result adequate as a nucle- 
ation agent. A concentration of 2% tetradecanol resulted sufficient 
to suppress subcooling both in tetradecane and microencapsu- 
lated tetradecane, suggesting that heterogeneous nucleation was 
more favourable in the presence of tetradecanol. Nevertheless the 
amount of nucleation agent should be the smallest possible in order 
to eliminate subcooling and reduce by the least the melting latent 
heat. Zhang et al. [19] also investigated crystallization from a DSC 
for the n-octadecane microencapsulated with l-octadecanol as a 
nucleation agent. The cooling curves obtained from the DSC for 
the microencapsulated n-octadecane were mainly affected by the 
average diameter of particles. The maximum degree of subcool- 
ing was 26°C with a cooling/heating rate of the DSC of 10°C/min. 
An addition of 10% in weight of l-octadecanol inside the micro- 
capsules as a nucleation agent decreased subcooling from 26°C to 
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Fig. 6. Effect of two different types of seeds on the nucleation process [42]. 


12°C with a DSC rate of 10°C/min. Yamagishi et al. [44] evaluated 
two n-tetradecane and n-dodecane microencapsulated slurries and 
observed that despite the fact that bulk PCM did not present 
subcooling, it was detected in the microcapsules. In the particle 
size range between 5 and 100 um, the crystallization temperature 
decreased along with particle size. |-Tetradecanol was added as 
nucleation agent and it cancelled subcooling. Griffiths and Eames 
[45] studied the behavior of a chilled ceiling through which amPCM 
slurry circulated. The DSC curves indicated that if the PCM was 
completely melted before being cooled, the mPCM slurry showed a 
subcooling of 5 °C. Nevertheless the experimental observations and 
the DSC indicated that if the mPCM slurry was cooled before com- 
plete melting of the PCM, the presence of solid material acted as a 
nucleus, and therefore the PCM solidified at higher temperatures. 
For this reason lower outlet temperatures than those of complete 
melting of the PCM were used. 


4.2. Stability 


4.2.1. Stability of PCM emulsions 

The stability question is always a concern during the develop- 
ment and production of an emulsion. There are five main instability 
problems [46]: 


e Creaming or sedimentation, as a result of gravity if there are den- 
sity differences between disperses and continuous phases. In an 
O/W emulsion (oil/water), creaming is the movement of the oil 
droplets due to gravity, causing a concentrated layer at the supe- 
rior part of the sample. Sedimentation is the opposite process. 


e Flocculation, the process of agglomerating the emulsion droplets. 

e Coalescence, process through which two or more droplets of the 
disperse phase melt and form a bigger droplet. 

e Ostwald ripening, as a result of solubility differences of the dis- 
perse phase contained in droplets of different sizes. 

e Phase inversion, process through which the continuous phase 
converts itself into the disperse phase and the disperse phase 
converts itself into the continuous phase. 


Creaming, sedimentation and flocculation are not accompanied 
by a change in size distribution of the disperse phase, while coa- 
lescence and Ostwald ripening cause an increase in the size of the 
droplet. Fig. 7 shows the instability processes mentioned. 

In a stable emulsion, the aforementioned instability processes 
never occur. However, the majority of emulsions inherently tend 
to separate, especially under complex operation conditions. Many 
studies analyze the stability of emulsions, which is of vital impor- 
tance in the future use of emulsions as thermal storage materials. 
It is also important that thermal properties remain invariable after 
thermal cycling. Huang et al. [10] observed that, after a one-month 
storage period of a RT 10 emulsion (paraffin from Rubitherm), 
emulsions with a 15-60% weight concentration presented cream- 
ing, while emulsions with a 60-75% concentration did not present 
instability. After one year, no changes in size distribution of par- 
ticles were detected and the emulsions maintained their heat 
capacity. Under mechanical-thermal loads, the droplet size of the 
paraffin emulsion suffered a slight increase in size. Besides, the 
nucleation agents of the paraffin droplets separated, showing a 
marked subcooling. The viscosity curves were similar before and 
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Fig. 7. Instability processes in emulsions [10]. 


after the tests and the heat capacity was maintained. Choi et al. [47] 
investigated the possibility of using an emulsifier to generate PCM 
droplets, specifically of hexadecane, which would not adhere to 
each other. The amount of emulsifier must be above a critical value 
in order to obtain emulsions with uniform and stable droplets. The 
emulsifier affected also the dynamics of the interface between an 
exchanger surface and the PCM, establishing whether hexadecane 
adhered to the surface. The project of Lorsch et al. [35] was carried 
out to determine the feasibility of applying a series of emulsions 
with paraffin mixtures hexadecane/tetradecane 70/30 in chillers, 
specifically to study the risk of clogging the pipes with these new 
heat transfer fluids. This agglomeration was more probable when 
the droplet size increased during the solidification processes. Pre- 
vious investigations have established that when the size of the 
droplet was below 0.1 mm, the pipes did not clog. This droplet size 
could be controlled through the amount of surfactant. Zou et al. 
[39] describe how agglomeration or precipitation processes can 
be avoided by: 1) reducing the diameter of the PCM droplets, as 
a very small droplet size causes a reduction in gravity force and 
the Brownian movement can be sufficient to compensate for grav- 
ity. The droplet size can be reduced through the optimization of 
emulsifying conditions, such as surfactant agents, their concentra- 
tion, agitation speed, duration of emulsification, etc. 2) searching 
for a reasonable concentration of PCM. Schalbart et al. [11] pre- 
pared several PCM nanoemulsions using different emulsification 
routes from low energy emulsification methods (phase inversion 
temperature method), determining the best emulsification routes 
to obtain satisfactory stability in the resultant emulsions. 


4.2.2. Stability of mPCM slurries 

In the case of mPCM slurries, problems such as sedimenta- 
tion, creaming or agglomeration can also occur, but the main 
problem is microcapsule rupture. Yamagishi et al. [44] studied 
the damage produced by the stresses caused by the pump or by 
agitation on microcapsules. They observed that agitation did not 


damage microcapsules and when suffering stresses caused by the 
pump, the rupture index decreased along with particle size. Regard- 
ing thermal cycling (contraction-expansion), it was observed that 
the 5 pum microcapsules did not rupture and the thermal prop- 
erties of cycled microcapsules remained invariable. Zhang et al. 
[20] studied the influential factors on thermal stability in rela- 
tion to the structure of n-octadecane microencapsulated into urea 
melamine formaldehyde capsules. A loss of liquid n-octadecane 
from the microcapsule was observed due to microcapsule rup- 
ture, as a consequence of a bad adjustment of core and capsule 
expansion at high temperatures. Zhang et al. [24] investigated the 
effects of several parameters of the microencapsulation process 
on the properties and stability of n-octadecane microencapsulated 
and nanocapsulated into melamine formaldehyde capsules. The 
thermal stability of the microcapsules and nanocapsules increased 
along with the increase in agitation speed and with the emulsi- 
fier content. Gschwander et al. [48] built an experimental setup 
to cycle different mPCM slurries. The mPCM slurries were pumped 
during several weeks in order to study stability. Conventional com- 
ponents such as expansion valves, heat exchangers, etc., worked 
well with mPCM slurries. The stability of microcapsules was neg- 
atively affected by high shearing speeds, which occurred mainly 
in the pump. After testing several pumps, the centrifugal pumps 
resulted in less damage or less destruction of the microcapsule 
shell. The good stability of capsules was a result of a smaller diam- 
eter and a thicker shell. This configuration resulted in a smaller 
fraction of paraffin and therefore a lower melting enthalpy. The 
microcapsules could be pumped during weeks, which resulted into 
a period equivalent to years in a real installation (800 cycles per 
day). Fan et al. [17] synthesized and studied in detail the thermal 
stability and permeability of microencapsulated n-octadecane and 
cyclohexane. The content of cyclohexane in the core has a remark- 
able effect on morphology, thermal stability and permeability. After 
thermal treatment of capsules, the volatile character of cyclohex- 
ane caused an expansible space to be formed, which improved the 
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thermal stability of microcapsules. These microcapsules presented 
a greater resistance to temperature and less mass loss. This weight 
loss was attributed mainly to the loss of n-octadecane in some 
broken capsules. Alvarado et al. [28] presented thermophysical 
properties of a microencapsulated tetradecane slurry. In concrete, 
with respect to microparticle stability, it was observed that the size 
of the microcapsule should be smaller than 20 um to avoid rup- 
ture. Ifa considerable amount of microcapsules were broken during 
durability tests, smaller capsules or capsules with a higher thick- 
ness/diameter ratio were tested. Microcapsules with a 2-10 um 
diameter showed less damage during circulation through a pro- 
gressive pump. Jin et al. [21] prepared and studied the behavior of 
paraffin microcapsules with a phase change temperature of 50.2 °C, 
specifically the influence of the weight percentage with respect to 
the PCM core on the stability during phase change. They observed 
that the microcapsules suffered coalescence when the content of 
the capsules was inferior to 2.1% in weight as a consequence of 
microcapsule rupture. When the content of the capsule was 16.7% 
in weight, the microcapsules suffered partial coalescence. When 
this percentage increased to 28%, the capsules remained intact after 
an elevated number of thermal cycles and presented a better ther- 
mal stability than bulk paraffin. Alkan et al. [22] characterized the 
behavior of docosane microencapsulated into PMMA (polymethyl 
methacrylate) capsules, which were synthesized by them. Degrada- 
tion of bulk docosane took place in a single stage at a temperature of 
approximately 120 °C, while microcapsules degraded in two stages: 
at 240°C and at 323°C. No significant changes were observed in 
the temperature or phase change enthalpies, after thermal cycles 
(1000, 3000, and 5000 cycles). Griffiths and Eames [45] pumped 
a mPCM slurry through a chilled ceiling. During the experiments, 
deposition of microcapsules was not observed in the pipe elbows 
and around valves. Degradation of the slurry was not observed 
either. It was observed that dry pumps were adequate to pump 
the slurry. 

From this review on the stability of PCM emulsions and mPCM 
slurries it can be concluded that, for the case of PCM emulsions 
and mPCM slurries, the creaming phenomenon and microcapsule 
rupture are the most frequent problems. In the case of emulsions, 
the creaming problem can be solved through the reduction of the 
PCM droplet size. In the case of the mPCM slurries, microcapsule 
rupture is solved through the reduction of particle size and/or 
increase of the relationship thickness capsule/thickness core for 
the PCM. 


4.3. Rheological behavior, viscosity and pressure drop 


The viscosity measurement and rheological characterization of 
these slurries systems, PCM emulsions and mPCM slurries, must 
be taken into account as although it seems homogeneous at a 
macroscopic level, the presence of other phases at a microscopic 
level modifies their behavior significantly. The concept of viscos- 
ity is based on the supposition that fluids are homogeneous and 
that, as a consequence, can be treated as continuous. Neverthe- 
less, disperse systems exhibit a behavior that diverges strongly 
from the behavior observed in homogeneous fluids. This is due 
to the presence of two or more phases, which generates one or 
more interfaces of different nature. As a consequence, the behavior 
of the flow can range from simply Newtonian in diluted particle 
systems with no interaction, to highly non-Newtonian behav- 
iors in concentrated or diluted systems where there is elevated 
interaction between particles. The main physical properties that 
influence rheological behavior are the content of the disperse 
phase; size, form and distribution of particle size and temperature. 
For the case of emulsions, two other factors must also be consid- 
ered: degree of droplet deformation and disperse phase viscosity 
[49]. 
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Fig. 8. Relationship between pumping work and heat transfer [50]. 


4.3.1. Rheological behavior of PCM emulsions 

Choi et al. [33] observed that in their 10% hexadecane emulsion, 
pressure drop decreased significantly starting from the melting 
point of PCM particles. The focus of the work of Royon et al. [5] 
was on the study of the thermorheological behavior of an emulsion 
with a 50% n-alkane concentration in weight. The results of the 
rheological experiments showed a pseudoplastic behavior and the 
apparent viscosity of the emulsion followed the power law. Vis- 
cosity decreased with the increase of temperature and the n and 
k indices were functions of temperature. The k index decreased 
with an increase in temperature, and the relationship between the 
k index and temperature was exponential. The n index showed 
a slight rise close to 9.5°C (phase change temperature), mainly 
caused by the phase change of the component in dispersion. Chen 
etal. [50] built an experimental installation to study the rheological 
characteristics of a tetradecane emulsion with a weight concentra- 
tion of 30%, prepared by the phase inversion temperature method. 
The emulsion could be considered a Newtonian fluid. The friction 
factor was adjusted to the classical function f=64/Re for laminar 
flow. Viscosity of the emulsion was 5.57 times that of water. The 
flow and pump work decreased considerably for the same amount 
of heat transfer when compared to water, due to phase change, 
as can be seen in Fig. 8. Pollerberg and Dötsch [51] presented the 
comparison of three slurries systems in distribution system appli- 
cations: ice slurries, tetradecane emulsion, and microencapsulated 
tetradecane slurry. One of the points studied was the rheological 
behavior of these slurries. The emulsion shows a Newtonian behav- 
ior in a weight concentration range of 10-20% and its viscosity was 
between 2 and 8 times that of water, depending on the tetradecane 
concentration. Both the slurry and the emulsion caused a higher 
pressure drop than water in laminar flow. Nevertheless, the pres- 
sure drop curves overlap in turbulent regime. All samples prepared 
by Huang et al. [10] presented a pseudoplastic behavior. The Kand n 
coefficient of the power law (t=K-") were practically constant with 
a weight fraction of 15-50%. A decrease in n and an increase in K 
were observed when concentration exceeded 50%. Besides, viscos- 
ity increased with the increment of solidified fraction. One reason 
for this behavior is that the solid droplets were not suffering sig- 
nificant deformation, resulting in a higher viscosity than when the 
droplets were melted. A working limit fraction was established at 
50%. Huang et al. [36] observed a pseudoplastic behavior in their 
PCM emulsions, in which viscosity decreased along with shear rate 
and then tended to remain constant when shear rate was higher 
than 30 s-!. Temperature influenced significantly, especially phase 
change. The reason for this phenomenon could be that the solid 
particles were not very deformed by shearing, causing a higher 
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viscosity than that of liquid droplets. Cho et al. [52] observed that 
the increase in pumping work in the case of their emulsion with 
a 25% volumetric concentration of PCM was small in the turbu- 
lent region, in comparison to water. In the case of the emulsion 
of Lorsch et al. [53], a fraction of volume up to 20% increased the 
pressure drop only by 3% compared to water. When the mixture 
was cooled to a temperature below PCM solidification, the pressure 
drop increased from 3 to 5%. Zou et al. [39] verified that for a given 
PCM concentration, there is a flow rate that causes the greatest 
savings in pumping consumptions. 


4.3.2. Rheological behavior of mPCM slurries 

Yamagishi et al. [44] studied the feasibility of n-tetradecane 
and n-dodecane microencapsulated slurries in low temperature 
applications. One of the studied aspects was the viscosity of 
these slurries. Specifically, viscosity was measured by a cylindri- 
cal Couette viscometer and it was observed that apparent viscosity 
depended on several parameters: PCM concentration, tempera- 
ture of slurry, and microcapsule size. When adding addictives 
such as surfactant agents, the slurry stopped behaving as a non- 
Newtonian fluid and started behaving as a Newtonian fluid. Yang 
et al. [14] prepared slurries with microencapsulated tetradecane 
from different materials and observed that the capsule material 
did not influence the viscosity of the sample. Rao et al. [54] cen- 
tered their study on the flow characteristics of microencapsulated 
n-octadecane slurries with concentrations from 5 to 20%, circu- 
lating through horizontal minichannels. The friction factor of the 
slurries under laminar regime increased along with the PCM con- 
centration. Compared to the friction fraction of water, a slight 
increase was observed in low-concentration (5%) slurries. Nev- 
ertheless, when concentration was 10% or higher, the increment 
in friction factor was more marked. The increment in PCM con- 
centration of the slurries tended to suppress the generation of 
turbulence in the flow. When concentration was at least 15%, no 
obvious transition was observed from laminar to turbulent flow as 
occurred in low-concentration PCM at Re= 2000. The pressure drop 
of the mPCM slurries flowing through the minichannels increased 
throughout the speed range when concentration increased. As pre- 
viously mentioned, Pollerberg and Dötsch [51] compared an ice 
slurry, a tetradecane emulsion and a microencapsulated tetrade- 
cane emulsion. The slurry presented the behavior of a pseudoplastic 
fluid, and viscosity was between 120 and 550 times that of water. 
Both the mPCM slurry and the PCM emulsion studied in this work 
presented a higher pressure drop than water in laminar regime. 
Nevertheless, the pressure drop curves overlapped in turbulent 
regime. In the microencapsulated tetradecane slurries developed 
by Alvarado et al. [28] with concentrations oscillating between 
5 and 17.7%, the relative viscosity seemed to be independent of 
temperature (in the phase change temperature range) for all con- 
centrations. The slurries behaved as a Newtonian fluid, until mass 
fractions of 17.7%. The results also indicated that the pressure 
drop increased slightly when PCM microcapsules were used, but 
did not affect significantly the pumping work. In some cases the 
pressure drop was lower than water, maybe due to a possible rup- 
ture of microcapsules and liberation of phase change material. The 
viscosity of the slurries prepared by Wang et al. [55] with concen- 
trations ranging from 5 to 27.6% of l-bromohexadecane adjusted 
to the values predicted by the Vand model. All slurries presented 
a Newtonian behavior. The pressure drop measurements showed 
a marked transition when the flow changed from laminar to tur- 
bulent flow. The friction factors in turbulent flow adjusted well 
to the classic Hagen Poseiuille model, while the friction factors 
in laminar flow were lower than those calculated with the Bla- 
sius equation, as observed in Fig. 9. In the experimental work of 
Heinz and Streicher [56], it was observed that with a concentra- 
tion up to 30%, the pressure drop was not considerably higher that 
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Fig. 9. Friction factors vs. Reynolds number [55]. 


water, for which this concentration resulted to be a good compro- 
mise between storage capacity and pressure drop. Chen et al. [30] 
observed that pumping power decreased considerably in compari- 
son to water due to the higher transported heat. The decrease could 
be of up to 67.5%. Zhang and Zhao [32] investigated the rheologi- 
cal properties of two mPCM slurries. Overall these mPCM slurries 
could be considered as Newtonian fluids when the shear rate is 
higher than 200s~! and PCM microcapsules concentration lower 
than 35%. The viscosity was higher for slurries with bigger PCM 
microcapsules. 

Table 4 shows the viscosity values of different PCM slurries stud- 
ied in literature. 

All consulted studies coincide when remarking that even with 
PCM concentrations close to 20-30%, the pressure drop of mPCM 
slurries and emulsions is slightly superior to water, and practically 
the same in turbulent regime. 


4.4, Thermal properties: Thermal conductivity 


One of the main disadvantages of thermal energy storage 
systems with PCMs is their low thermal conductivity. This circum- 
stance results in slow charging and discharging. In the bibliography, 
there are numerous studies aimed at the improvement of the ther- 
mal conductivity of the PCM, either by embedding structures of 
materials with high thermal conductivity, or by using finned heat 
exchangers or encapsulating the PCM in containers with a high 
surface/volume ratio. This is the reason why PCM microcapsules 
are interesting. Due to the microscopic size of the PCM microcap- 
sules or droplets, the PCM slurry can be treated as a homogeneous 
material. This assumption implies in that the temperature gradients 
inside the solid are negligible. This is accomplished if the convective 
thermal resistance inside the microcapsules is low in comparison 
to the convective thermal resistance between the microcapsule and 
its surroundings. The fulfilment of this condition can be analyzed 
through the Biot number, which should be under 0.1. If the Biot 
number is calculated for a PCM microcapsule, a value much lower 
than 0.1 will be obtained. This means that there is no conduction 
inside the microcapsule and therefore temperature gradients can 
be neglected. When the PCM is microencapsulated, low thermal 
conductivity ceases to be a problem [58]. 

So, PCM slurries in water can improve heat transfer as a conse- 
quence of the relationship area/volume of droplets in the case of 
emulsions and of microcapsules in the case of slurries, in compari- 
son to systems in which the PCM is macroencapsulated. Besides, the 
fact of dispersing phase change particles into a fluid can improve 
heat transfer through convection with respect to water. These slur- 
ries can serve either as thermal storage materials or heat transfer 
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Table 4 
Viscosity values for different PCM emulsions and mPCM slurries studied in literature. 
PCM Weight concentration Viscosity Ref. 
Mixture of n-alkanes 50% 25 mPas (7°C, stress = 3 Pa) [5] 
10% 1.25 mPas (35°C 
Tetradecane 20% 2.33 mPas (38°C [14] 
PCM emulsion 30% 2.95 mPas (38°C 
5% 2.05 mPas (10°C 
1.54 mPas (20°C 
l-Bromohexadecane 40% 2.23 mPas (10°C [57] 
: . 1.81 mPas (20°C 
Dispersion system 15.80% 4.18 mPas (10°C 
3.32 mPas (20°C 


Newtonian behavior 
starting at 600 1/s 


5% 0.98 mPas (25°C 
2.74mPas (5°C 
10% 1.04mPas (25°C 
1.85 mPas (5°C 
Tetradecane 15% 1.54mPas (25°C [14] 
mPCM slurry 1.57 mPas (5°C 
20% 1.94 mPas (25°C 
1.39 mPas (5°C 
25% 3.31 mPas (25°C 
5.23 mPas (5°C 
30% 4.21 mPas (25°C 
5.53 mPas (5°C 
35% 4.11 mPas (25°C 
5.84 mPas (5°C 
40% 12.58 mPas (25°C) 
20.87 mPas (5°C 
5% 1.27 mPas 
n-Octadecane 10% 2.3 mPas [29] 
20% 4.9 mPas 


Water 


1.307 mPas (10°C) 
1.002 mPas (20°C) 
0.798 mPas (30°C) 


fluids. The thermal properties of these slurries are different from 
those of PCM and the fluid in question, which are essential to eval- 
uate the fluid and the heat transfer characteristics of a system with 
these slurries. The thermal properties to be discussed are thermal 
conductivity and convection heat transfer coefficient. The analysis 
of the different studies regarding the convection heat transfer coef- 
ficient is presented in a separate section, due to their extension and 
importance within this review. 

The thermal conductivity of diluted dispersions, Ag, can be eval- 
uated from Maxwell’s relationship: 


Aq 24+ (Ap/kp)+2-€-(Ap/(Ap - 1)) 
Ap 2+(àp/kf)— c- (Ap/(Ap - 1) 


where Ap is the thermal conductivity of the disperse phase, A, is 
the thermal conductivity of the continuous phase and c is the vol- 
umetric concentration of the disperse phase. As a consequence of 
the interactions particle/fluid, the effective thermal conductivity is 
greater than that predicted by this equation, and can be obtained 
from the following correlation: 


(1) 


— =1+B-c-Pem (2) 


where Pe, is the Peclet number for the particle. The values of B and m 
depend on the Peclet number of the particle. In general, the thermal 
conductivity of phase change particles is not very elevated and must 
be improved [59]. Xuan et al. [60] developed a new type of heat 
transfer fluid, consisting of a magnetic slurry with microencapsu- 
lated PCM that incorporated the advantages of microencapsulated 
PCM and of a magnetic fluid. Iron nanoparticles were added to the 
melamine-urea-formaldehyde capsule surrounding the paraffin. 


Thermal conductivity was measured through the hot wire method. 
The measured data pointed out that among all particle components, 
the content of iron nanoparticles was the dominating factor on ther- 
mal conductivity. The thermal conductivity ratio increased linearly 
with the increase in concentration of iron nanoparticles inside the 
capsules. Similarly, Ho and Gao [61] prepared a n-octadecane emul- 
sion with alumina particles. Compared bulk paraffin, the emulsion 
presented a higher thermal conductivity, depending on the concen- 
tration of the alumina particles. This relationship was not linear. 
Improvements between 2 and 6% were obtained with weight con- 
centrations of alumina between 5 and 10%. Ho et al. [62] managed 
to improve the thermal conductivity of PCM emulsions by adding 
alumina nanoparticles, to values even better than those of water. 
Table 5 shows some thermal conductivity values for studied slurries 
in literature. 


5. Heat transfer 


The main objective of the publication by Kasza and Chen [63] 
was to point out the benefits of the use of PCM slurries in water, 
such as the improvement in heat transfer and increase in storage 
efficiency. Some of the specific benefits mentioned were: 1) reduc- 
tion in the temperature difference between source and drain, 2) 
increase of the heat capacity of the fluid, as a consequence of the 
PCM dispersion, giving place to a lower mass flow and therefore a 
lower pumping consumption, 3) dynamic use of the PCM. In a con- 
ventional system, heat exchange between PCM (static use) and a 
separated heat transfer fluid is needed to transport heat or cooling. 
Nevertheless with PCM slurries, thermal storage and the heat trans- 
fer fluid are integrated into the PCM slurry. The energy losses of the 
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Table 5 
Thermal conductivity values for mPCM slurries studied in literature. 


PCM Weight concentration Thermal conductivity (W/(m K)) Ref. 
5% 0.571 (a 20°C) 

n-Octadecane 10% 0.541 (a 20°C) [29] Calculated with Eq. (1) 
20% 0.483 (at 20°C) 
5% 0.568 (at 20°C) 

mPCM slurry 10% 0.539 (at 20°C) 

1-Bromohexadecane 15.80% 0.506 (at 20°C) [55] Calculated with Eq. (1) 
20.40% 0.48 (at 20°C) 
27.60% 0.446 (at 20°C) 


n-Eicosane 


Agua 


5% (with 10% alumina nanoparticles) 
10% (with 10% alumina nanoparticles) 


0.66 (at 30°C) 
0.62 (at 30°C) [62] 
0.60 (a 20°C) 

( 


0.61 (a 30°C) 


heat exchange will be annulled, as the necessity for a secondary 
heat transfer fluid is eliminated. 

The possible improvement mechanisms for heat transfer are 
enumerated. Improvement in heat transfer occurs in slurries, with 
or without phase change. This improvement is substantially greater 
when considering PCM slurries. Several mechanisms responsible 
for this improvement were examined, including particle rotation 
and migration in the flow as well as the role of phase change. 

The following sections show a review of the studies carried out 
up to date on the heat transfer phenomenon (mainly forced con- 
vection inside a tube) and were classified in accordance to the 
type of fluid (mPCM slurry or PCM emulsion), to the character 
(experimental or numerical), and to flow regime (laminar or turbu- 
lent). Boundary conditions and geometry were also accounted for. 
The following compilation aims at completing the subsection on 
numerical models of heat transfer phenomenon in mPCM slurries 
of the review of Dutil et al. [64], on PCM and mathematical models. 


5.1. Heat transfer by forced convection in PCM emulsions 


5.1.1. Laminar flow. Experimental studies 

Roy and Avanic [65] presented the results of an experimen- 
tal study on heat transfer by convection in a circular tube, 
through which circulated a n-octadecane emulsion. The experi- 
mental results showed that the characteristics of heat transfer by 
convection in laminar flow for the emulsion were similar to those 
of microencapsulated PCM slurry. This means that the walls of 
the microcapsules did not have a significant effect on the heat 
transfer process. The results also showed that heat transfer as a 
whole did not change significantly even with very high concen- 
trations. The experiments also verified that the Reynolds number 
was not an independent parameter for the heat transfer process. 
From an application viewpoint, subcooling and the melting tem- 
perature range must be considered in future works. The type of 
slurry (PCM emulsion or mPCM slurry) must be selected consider- 
ing fluid-dynamics criteria and not heat transfer criteria. 


5.1.2. Turbulent flow. Experimental studies 

Choi et al. [33] investigated the increase of the convective heat 
transfer coefficient in the case of a hexadecane emulsion in an 
experimental installation. The local convection heat transfer coeffi- 
cient improved with respect to water. The study proposed a melting 
model with three regions: a first region where the average emul- 
sion temperature is below melting temperature, a second region 
where phase change takes place, and a third region where the PCM 
is completely melted. The local convective heat transfer coefficient 
increased in region 1, decreased abruptly in region 2, and increased 
slightly in region 3. The authors use the hypothesis that region 1 
presented a very thin layer of melted particles, which increased 


rapidly in region 2 preventing the transportation of solid parti- 
cles to the wall, and in region 3 there were no solid particles. As 
average temperature increased in the direction of the flow, the vis- 
cosity decreased and the Reynolds number increased, giving place 
to a slight increment in the heat transfer coefficient. It is also spec- 
ulated that when the heat flux is very high, the layer of melted 
particles close to the wall becomes thicker, preventing the radial 
migration of solid particles from the center to the wall, resulting in 
this way ina lower convection coefficient. Choi and Cho [34] studied 
experimentally the effect of the ratio height/width of rectangular 
channels through which circulated a paraffin emulsion (C22H46), of 
its concentration, of the heat flux and of the Reynolds number, on 
the heat transfer phenomenon. The local heat transfer coefficients 
were higher for emulsions than for water as a consequence of: 1) 
rupture of the laminar sublayer by the paraffin particles, 2) direct 
heat transfer due to the impact particle-wall, and 3) latent heat 
of paraffin. The influence of the ratio height/width of the channel 
on the heat transfer coefficients both for water and emulsion were 
higher with higher heat fluxes. 


5.2. Heat transfer by forced convection in mPCM slurries 


5.2.1. Laminar flow. Numerical studies 

Charunyakorn et al. [59] modeled the heat transfer phenomenon 
for mPCM slurries circulating through a circular tube under a 
constant heat flux. The mathematical model was solved using 
an implicit finite difference method. Heat generation (or absorp- 
tion) was included in the energy equation as a heat source due 
to phase change of the particles. The Stefan number and the PCM 
concentration were the dominating parameters. The effect of mod- 
ifying the Peclet number was not significant. The effect of the ratio 
of duct/particle radiuses was weak but perceptible in the range 
50-200 and insignificant for a radius range 200-400. The results 
showed that the slurry could improve the heat transfer coefficient. 
This improvement could be used as a higher heat transfer rate, a 
reduction in mass flow, or a more effective control of the wall tem- 
perature. Zhang and Faghri [66] presented a numerical solution for 
the problem of convection heat transfer in mPCM slurry in a cir- 
cular tube with constant heat flow under laminar flow conditions. 
The authors pointed out that one of the most important reasons for 
the differences in heat transfer between previous experimental and 
theoretical studies could be the non-consideration of subcooling, 
the microcapsule coating the PCM, and the melting temperature 
range. Their numerical model considered these factors and in this 
way, the difference between numerical and experimental results 
was decreased. The results also showed that heat transfer by forced 
convection was reduced with an increase in width of the phase 
change temperature range. Roy and Avanic [67] introduced the spe- 
cific heat as a function of temperature in the energy equation when 
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studying the heat transfer phenomenon in mPCM slurries, under 
laminar regime conditions and constant heat flux. The effects of 
phase change were incorporated directly in the energy equation, 
assuming the heat capacity of the phase change material as a func- 
tion of temperature. The numerical study showed that the Stefan 
number was the only parameter with a significant impact on heat 
transfer for the typical mPCM slurries and for typical operation 
conditions. The subcooling effects could be important at low heat 
fluxes or when the input temperature was much lower than melt- 
ing temperature. The melting temperature range had a significant 
impact when it was wide. A correlation was developed to pre- 
dict the increment in wall temperature in function of the distance 
along the length of the tube. Hu and Zhang [68] studied numer- 
ically the convection heat transfer of mPCM slurries circulating 
through a circular tube under constant heat flux. It was concluded 
that the exact nature of the phase change process affected the 
improvement degree of heat transfer through convection. The Ste- 
fan number and the PCM concentration resulted to be the most 
influential parameters on the improvement of heat transfer in 
mPCM slurries. Nevertheless, subcooling, phase change tempera- 
ture range and particle diameter also influenced the improvement 
of heat transfer. Heat transfer improved when subcooling and phase 
change temperature range decreased, and when particle diameter 
increased. The improvement would be more notable in a tube with 
smaller diameter and greater in the completely developed region 
than in the entry region. Lu and Bai [69] presented a numerical 
model to analyze heat transfer by convection in mPCM slurries 
under laminar flow. When the microcapsule radius increased, the 
heat transfer efficiency increased. The increase in the volumetric 
fraction of the PCM particles could result in an improvement of 
heat transfer for a constant heat flux. When the Stefan numbers 
increased, the heat transfer efficiency was lower. Heat transfer also 
improved when the Reynolds number was increased. Xing et al. 
[70] carried out a parametric study from a numerical model for 
mPCM slurries in microchannels. From this study it was concluded 
that for a given Reynolds number and PCM concentration, there 
was an optimum heat flux through the wall that resulted in a max- 
imum effectivity factor (ratio between the heat absorbed by the 
mPCM slurry and the heat absorbed by water). Also for a given 
Reynolds number and PCM concentration, there was a maximum 
index for the optimal heat flux; where this index compared the 
ratio between the absorbed heat and pumping work of the slurry 
against water. When the Reynolds number increased, the maximal 
effectivity factor decreased. In general, the increase in Reynolds 
resulted in a decrease of the ratio heat absorbed/pumping power. 
The model also determined that the results were very sensible to 
the inlet temperature of the fluid in the test tube, which could 
explain the discrepancies with the experimental results of liter- 
ature. The convection heat transfer coefficient reached a peak at 
the moment of phase change. In the region before and after melt- 
ing, the heat transfer coefficients were lower than water. Fig. 10 
shows the local heat transfer coefficient for the slurry for different 
Reynolds numbers for a given heat flux. Zhao et al. [71] presented 
the results of simulations to evaluate the convective heat trans- 
fer for mPCM slurries, circulating through a circular tube with 
the boundary condition of constant wall temperature. The most 
influential parameters on the improvement of heat transfer were 
the Stefan number, the volumetric concentration, and the ratio 
particle/tube radiuses. Heat transfer improved when the Stefan 
number decreased, when volumetric concentration increased or 
when the ratio particle/tube radiuses increased. The influence of 
particle diameter on the improvement of heat transfer originated 
from the microconvection caused by the microcapsules in sus- 
pension, instead of the high heat capacity due to phase change. 
The Reynolds number was the most dominating parameter on the 
length of the phase change region. The authors remark that it is 
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Fig. 10. Heat transfer coefficients for different Reynolds numbers [70]. 


difficult to develop a simple and general correlation due to the 
complexity of the combined effects. Sabbah et al. [72] used Fluent 
to solve the heat transfer phenomena in rectangular microchan- 
nels (width 100 um, height 500 um and length 10 mm), through 
which circulated a mPCM slurry in laminar flow under constant 
heat flux conditions. The use of low concentrations of PCM gave 
place to a better cooling than the case of water for high heat flux. 
The improvement index was greater for lower concentrations due 
to the elevated increase in the viscosity of slurry along with concen- 
tration. The use of the mPCM slurry against water only reduced the 
temperature of the heat dissipated. Increasing the velocity of the 
slurry increased the heat transfer coefficient, always when the PCM 
melted inside the channel. In order to achieve a high improvement 
index for a given heat flux, the cooling system must be designed 
so that the PCM particles start melting at the entrance of the chan- 
nel and, when leaving, are completely melted. The use of the slurry 
properties dependent on temperature was crucial. Zeng et al. [57] 
adopted both experimental and numerical methods in their study, 
to analyze the convective heat transfer coefficient in microencapsu- 
lated l-bromohexadecane slurries in a circular tube under constant 
heat flux conditions. This work considered the heat capacity as a 
sine function of temperature. The investigation focused on the com- 
pletely developed thermal region to avoid the influence of thermal 
development. Three types of fluid were studied: pure water, water 
with microcapsules without phase change, and mPCM slurry, to 
separate the influence of phase change and microcapsules. The Nus- 
selt number and the dimensionless wall temperature were utilized 
to study heat transfer. The addition of microcapsules without phase 
change improved effective thermal conductivity. This improve- 
ment of the thermal conductivity resulted in an improvement of 
the local Nusselt number. With the PCM microcapsules, the Nus- 
selt number tended to draw a wave in the phase change region, 
as can be observed in Fig. 11. However, the average Nusselt num- 
ber did not change very much. During phase change, the authors 
attribute the decrease in wall temperature mainly to the contri- 
bution of the variable Cp, and not to the improvement in the heat 
transfer index between flow and wall. The amplitude of fluctuation 
of the Nusselt number increased with the decrease in phase change 
temperature range. This fluctuation also increased when the Ste- 
fan number decreased during phase change. The Nusselt number 
increased when the PCM concentration also increased, along with 
increases of the Reynolds number and the microcapsule diame- 
ter. Kuravi et al. [73] developed a 3D model to simulate the heat 
transfer phenomenon in a nanoencapsulated octadecane slurry in a 
poly alpha olefin-based oil (dielectric oil use in cooling) circulating 
through microchannels (width 200 um, height 50 um and length 
250 um) in laminar regime and under boundary conditions of 
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Fig. 11. Calculated Nu and temperature profile of the wall and PCM suspension [57]. 


constant heat flux. A sine function was used to model phase change. 
The results showed an increase of the Nusselt number with an 
increase of PCM concentration, which decreased with the increase 
in heat flux. For the considered mass flows, an increment in heat 
flux resulted in an improvement of the heat transfer phenomenon. 
The authors comment that, based on their numerical results, if 
the melting temperature range is narrow, the temperature at the 
microchannel entry should be close to the peak temperature of the 
PCM melting curve. These results depended on the dimensions of 
the channel. 


5.2.2. Laminar flow. Experimental studies 

Goel et al. [74] carried out an experimental study utilizing 
a microencapsulated n-eicosane slurry to evaluate heat transfer 
under laminar flow conditions in a circular tube with constant heat 
flux. It was concluded that the Stefan number was the dominating 
parameter in heat transfer. The volumetric concentration did not 
have a significant direct effect on heat transfer. Nevertheless, con- 
centration is included in the definition of the Stefan number and 
therefore has an indirect effect. A higher ratio between the diame- 
ters microcapsule/tube improved heat transfer. It was also proved 
that the homogeneity degree of the slurry could be neglected, but 
required a deeper experimental investigation. The results adjusted 
well qualitatively with the work of Charunyakorn et al. [59], but 


the quantitative adjust was not good. The results showed that a 
decrease of 50% could be achieved in the temperature wall using 
microcapsules. The best results were achieved utilizing lower Ste- 
fan numbers (Ste < 1), i.e., either with a PCM with a high latent heat, 
or increasing concentration, or reducing the diameter of the tube. 
The authors also remarked that another important parameter was 
the inlet temperature or initial subcooling. For inlet temperatures 
lower than melting temperature, the reduction in temperature due 
to energy storage was eclipsed by the sensible heat process until 
melting temperature is reached. Rao et al. [29] compared experi- 
mentally the results obtained using water and microencapsulated 
n-octadecane slurries with weight concentration from 5 to 20%, 
under the same thermal conditions and mass flows, when cir- 
culating through rectangular microchannels (hydraulic diameter 
2.71 um) under laminar regime. The mPCM slurry showed a better 
cooling behavior compared to water for low mass flows. The reduc- 
tion in the wall temperature was higher when the concentration 
of microencapsulated PCM increased, and the heat transfer coef- 
ficient and Nusselt number increased with concentration. When 
the mass flow was increased, the cooling of the mPCM slurry was 
less effective. For each slurry, the heat transfer coefficient and 
the Nusselt number increased along with the increase of mass 
flow, but this increase was less noticeable when the PCM con- 
centration increased. Heat transfer at higher concentrations was 
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less effective in high mass flows due to a lower residence time 
in the channels, lower thermal conductivity and less active move- 
ment of the slurry particles. The experiments of Wang et al. [55] 
to study heat transfer in a circular tube covered both the case of 
laminar flow and slightly turbulent flow for a microencapsulated 
l-bromohexadecane slurry under constant heat flux. The local heat 
transfer coefficients obtained for the mPCM slurry were higher than 
water. The authors attribute this improvement to the latent heat 
associated with phase change. In the case of laminar flow, these 
local heat transfer coefficients increased with the fraction of PCM 
in the fluid. In the case of turbulent flow, the local heat transfer 
coefficients were significantly influenced by the heat flux through 
the test section, and the maximum values of the heat transfer coef- 
ficients occurred before for a higher heat flux. Heat transfer was 
considerably influenced by the turbulence degree of the flow. A 
correlation for laminar regime was proposed, which predicted sat- 
isfactorily the local heat transfer coefficient with a 15% error. 


5.2.3. Turbulent flow. Numerical studies 

Roy and Avanic [75] modeled heat transfer under turbulent flow 
conditions in a circular tube with constant heat flux, as previous 
publications only presented numerical models considering laminar 
flow conditions while most engineering applications involve turbu- 
lent flows. In the proposed method, the effects of phase change are 
directly incorporated in the energy equation, assuming heat capac- 
ity as a function of temperature. The most influential parameter in 
heat transfer was found to be the Stefan number. A considerable 
reduction in wall temperature was observed for low and moderate 
values of the Stefan number. The melting temperature range and 
the subcooling degree were also influential parameters. Royon and 
Guiffant [6] described a model to study heat transfer in a circular 
tube for a slurry with millimetric PCM particles, under constant 
wall temperature conditions. The model is useful to optimize the 
size of a heat exchanger and allows for the definition of the min- 
imum length of the heat exchanger to store a specific amount of 
latent energy from the particles. An increase in the mass fraction 
of particles increased the minimum length from a linear relation- 
ship. The effect of the Reynolds number on the minimum length of 
the heat exchanger was small. A correlation was proposed from the 
numerical results to determine the characteristics of the exchanger 
for a PCM slurry. 


5.2.4. Turbulent flow. Experimental studies 

Alvarado et al. [28] also presented experimental results on the 
convective coefficient in a circular tube under constant heat flux 
in turbulent regime. The heat transfer experiments showed that 
the improvement in heat capacity was considerable even at low 
mass fractions. The heat transfer coefficient, however, was lower 
for the case of slurry than for the case of water at the same velocity. 
An improvement in the thermal conductivity of the slurry could 
improve heat transfer. The heat transfer coefficient increased con- 
siderably during the phase change process. Fig. 12 shows the heat 
transfer coefficients for a slurry with 7% PCM weight concentration. 

Table 6 shows a summary of the studies carried out on the heat 
transfer phenomena by internal convection in these slurry systems, 
according to the slurry flow regime, the type of study, geometry and 
boundary conditions. 


5.3. Heat transfer through natural convection in mPCM slurries 


Heinz and Streicher [56] studied experimentally a 200-liter tank 
with a mPCM slurry developed by BASF with a melting tempera- 
ture of 60°C, with a typical spiral exchanger, inside which water 
circulated. As the limiting factor for heat transfer is the natural 
convection of the surface of the exchanger to the storage fluid, it 
resulted interesting to study the heat transfer coefficient by nat- 
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Fig. 12. Heat transfer coefficient for a suspension with 7% PCM [28]. 


ural convection. The heat transfer coefficients decreased as phase 
change process takes place, due to the reduction in the difference 
of temperatures between the exchanger and the storage fluid. Due 
to higher viscosities, the heat transfer coefficient through natu- 
ral convection decreased when the PCM concentration increased 
in water. Even with the lowest PCM concentration, 20%, the val- 
ues of the heat transfer coefficients for natural convection were 
lower than for water. A very similar study was carried out by Dia- 
conu et al. [78]. An experimental study was carried out on the heat 
transfer phenomenon of natural convection in a tank filled with a 
microencapsulated RT6 slurry (organic PCM from the manufacturer 
Rubitherm), where the heat exchanger through which water circu- 
lated consisted of a helicoidal copper tube, used as storage tank for 
solar air conditioning applications. During phase change of PCM, it 
was observed that the natural convective coefficient could be up 
to five times that of water, depending on the temperature condi- 
tions. The reason for this improvement could not be justified, as 
the phase change temperature range overlapped the temperature 
range in which water presented a drop in the natural convection 
coefficient. Huang et al. [79] also studied a heat storage system with 
a helical coil heat exchanger using a mPCM slurry with a phase 
change temperature 65°C (produced by BASF) with a 25, 35 and 
50% of volumetric concentration of PCM. The mPCM slurry was 
used as thermal energy storage material and water as secondary 
heat transfer fluid. The results showed that the mPCM slurry with 
a volumetric concentration of 50% was not suitable, since the low 
thermal conductivity and high viscosity reduced the heat transfer 
from the heat exchanger to the storage media. 


5.4. Other studies 


Heinz and Streicher [56], besides studying the phenomenon of 
convection in a tank filled with a slurry, also studied the heat trans- 
fer water-mPCM slurry in a plate heat exchanger. It was observed 
that the convective heat transfer coefficient decreased with the 
concentration of PCM, due to a greater viscosity and lower ther- 
mal conductivity. With a 20% concentration it decreased 30% in 
comparison to water, and with a 40% concentration, a decrease 
of 40% was observed. In spite of the lower overall heat transfer 
coefficient, the exchanged power of the plate heat exchanger is 
approximately the same as that of water. This is due to the higher 
difference in average temperatures, as result of the higher heat 
capacity of the slurry. The majority of the works, both experimental 


Table 6 
Compilation of studies carried out on the heat transfer phenomenon in PCM emulsions and mPCM slurries. 
Flow regime Type of study Boundary condition Geometry Simulation Experimental validation? Ref. 
tool/Numerical 
formulation 
Laminar Constant heat flux Circular tube - - [65 
PCM emulsions Experimental Constant heat flux Circular tube - - [33 
Turbulent Constant heat flux Rectangular - - [34 
channels 
Constant heat flux Circular tube Finite differences. 2D Ahuja [76]. Experimental validation [59 
without phase change. Differences not 
quantified. 
Constant heat flux Circular tube Finite differences. 2D Goel et al. [74]. Differences of 34%. [66 
. Constant heat flux Circular tube Fortran 90. Finite Roy and Avanic [67]. Very small [67 
Numerical differences. 1D differences. 
i Constant heat flux Circular tube Finite differences. 2D Goel et al. [74]. Differences under 6%. [68 
Laminar Constant heat flux Circular tube Finite differences. 2D Goel et al. [74]. Numerical results [69 
adjusted well to experimental results. 
mPCM slurries Constant heat flux Circular tube Finite differences. 2D Goel et al. [74]. Good adjust between [70 
experimental and numerical results, 
and sensible to entry temperature. 
Constant wall temp Circular tube Finite differences. 2D No. All available experimental results [71 
were limited to constant heat flux 
conditions. 
Constant heat flux Rectangular Fluent 6.2. 3D Goel et al. [74]. Adapted the geometry [72 
channels of their model. The numerical results 
adapted well to the experimental 
results. 
Constant heat flux Circular tube Finite differences. 2D Zeng et al. [57]. Differences under 9.4%. [57 
Constant heat flux Rectangular Comsol. Finite Goel et al. [74]. Adapted the geometry [73 
channels elements. 3D of their model. The numerical results 
adapted well to the experimental 
results. 
Constant heat flux Circular tube - - [74 
Experimental Constant heat flux Rectangular - - [29 
channels 
Constant heat flux Circular tube - - [55 
. Constant heat flux Circular tube Finite differences. 2D Choi [77]. Differences around 10-20%. [75 
Numerical Constant wall temp Circular tube Mathematica. Finite Validation with analytical solution [6] 
Turbulent differences. 1D without phase change. 
Experimental Constant heat flux Circular tube - - [28 
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and numerical, have assumed that the dispersed PCM and water 
are in thermal equilibrium. Nevertheless, Diaconu [80] analyzed 
through simulation the heat transfer phenomenon between the 
PCM microcapsules and water. The results considered a convective 
coefficient between microcapsules and water, and it was observed 
that the temperatures of water and PCM microcapsules were very 
close, which was attributed to the great exchange area between 
the PCM microcapsules and water. The greater differences between 
the temperatures of water and PCM microcapsules occurred dur- 
ing melting/solidification. Inaba et al. [81] developed a numerical 
two-dimensional model to study natural convection in a rectan- 
gular building with a non-Newtonian slurry of mPCM. The heat 
transfer coefficients were higher when compared to slurries with- 
out phase change. When the PCM concentration increased within 
the range 20-40%, the convection coefficient decreased, while for 
the range 10-20%, the convective coefficient improved. 


6. Applications 


For the moment, the main application present in literature is the 
utilization of these PCM emulsions and mPCM slurries as thermal 
storage materials and heat transfer fluids in chilled ceilings. Wang 
and Niu [82] presented the results of a mathematical simulation of 
acombined system of chilled ceiling and storage tank with amPCM 
slurry, in addition to an air treatment unit for the ventilation neces- 
sities, in a room with the climatology of Hong Kong. The slurry was 
cooled and stored in the tank during the night, which resulted in 
electricity peak shaving, taking advantage of the nocturnal tariff 
and of a higher COP of the machine due to operation during lower 
environmental temperatures. During working hours, the mPCM 
slurry flowed from the tank to the chilled ceiling, melting the 
PCM and releasing the latent heat. The combination of the chilled 
ceiling plus storage tank against a conventional water system 
achieved peak shaving, and therefore a smaller cooling unit/chiller 
could be sufficient. Three systems were compared regarding energy 


Table 7 


consumption: 1) Conventional system without storage, 2) Noctur- 
nal storage with mPCM slurry, 3) Storage with ice, where the loads 
of the air treatment unit and chilled ceiling were supplied by the 
tank. The consumptions were practically the same for the mPCM 
slurry and water. Nevertheless, it must be taken into account that 
calculations were carried out using the same COP for the case of 
the tank with water and mPCM slurry, when in reality the COP for 
the case of the tank with mPCM slurry is higher due to operation 
at lower environment temperatures (charging during the night). 
The electrical power will be lower. For ratios up to 2.4 (diurnal 
tariff/nocturnal tariff), the tank with mPCM slurrry presents higher 
economic benefits. Griffiths and Eames [45] studied experimentally 
the pumping of a mPCM slurry from BASF manufacturer through a 
chilled ceiling in a room. The room was tested during four months 
with a 40% PCM concentration. When water was pumped through 
the chilled ceiling, a mass flow of 0.71/s was required for an inlet 
temperature of 16°C and outlet temperature of 18 °C, maintaining 
the room at 19°C. When water was substituted by the mPCM slurry, 
the slurry was capable of maintaining a temperature of 20-21 °C 
with a mass flow of 0.251/s. This means that the ceiling required 
a lower mass flow (pumping savings were not quantified), could 
absorb energy at a constant temperature, avoiding increments in 
the panel surface temperature when internal gains increased. 
Another well known application, similar to the previously 
described, was carried out at the Narita Airport in Tokio by Shibu- 
tani [83]. The issue in the installation of the Narita Airport in Tokio 
was the change of refrigerants due to environmental reasons. When 
R11 and R22 were substituted by R134a and R123 without chang- 
ing the chiller unit, this resulted in lower cooling power and the 
chiller was non-capable to absorb the demand peaks at specific 
times of the day. This problem was solved through the installation 
ofa tank filled with amPCM slurry custom-developed by Mitsubishi 
Heavy Industries. The characteristic temperatures on the demand 
side were a supply temperature of 5°C and a return flow temper- 
ature of 12°C. A mPCM slurry was selected with a phase change 


Objective magnitudes and influential parameters at the time of selection of a PCM emulsion or mPCM slurry as heat transfer fluid or thermal storage material. 


Influential factors or parameters Objective magnitudes 


Influence when the factor increases 


Positive influence 


Negative influence 


Rupture of microcapsules 


Rupture pressure of microcapsules decreases, 
higher number of ruptured capsules. 


Particle diameter 7 
Subcooling 


Greater probability of existing nucleation 


agents, and therefore lower subcooling. 


Apparent hysteresis 


Possible non-equilibrium between PCM and water 
temperatures, possibility of hysteresis 


Heat transfer 


Improvement in convection coefficient. 


Stability of emulsions 


Heat capacity 


Creaming speed increases 


Increase in heat capacity, increase in 
transported heat. 


PCM concentration 
Pressure drop 


Increase of viscosity, increase of pressure loss and 
pumping work. Up to PCM concentrations of 
15-20% the increase is slightly superior to water. 


Decrease in Stefan number and therefore 


Heat transfer 


improvement of convection coefficient. 


Increase in viscosity, decrease in turbulence 
degree, and therefore worsening or convection 
coefficient. 


Decrease of thermal conductivity, occasioning 
deterioration in heat transfer. 


Operation temperature range 


Heat transfer 


The operation temperature range must fit 


with the phase change temperature range, 


and be the narrowest possible. 
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temperature range between 5 and 8 °C. The demand peaks occurred 
between 8:00 and 22:00, and therefore the cooling produced dur- 
ing the night by the chiller unit could be stored and reduce the 
demand peaks during the day. The slurry presented a storage den- 
sity of 67 MJ/m3, lower in comparison to an ice tank, 167 MJ/m3, 
but higher in comparison to water, 21 MJ/m3. Both the COP of the 
system and the operational costs for water and mPCM slurry were 
similar and lower than the ice tank. 

Pollerberg and Détsch [51] proposed an emulsion with a 20% 
weight concentration of tetradecane (melting temperature =5 °C) 
for cooling supply networks. In this way, the required volumetric 
flow was low, allowing for the reduction of the pumping power 
and pipe dimensions, with lower operation and investment costs. 
The results compared three scenarios for a fictitious cooling supply 
network: tetradecane emulsion, ice slurry and water. The calcula- 
tions considered the generation and distribution costs of cooling. A 
power of 500 kW was assumed, with a distribution length of 600 m 
and 1000 operation hours throughout the year. 

As a consequence of the higher heat capacity of the ice slurry, 
these presented lower distribution costs. The distributions costs 
for the emulsion were lower than those of water. The calculation 
of distribution costs takes into account the rheological behavior at 
the time of calculating the pressure drop. Due to the fact that the 
temperature of the evaporator must be lower in order to produce 
ice slurries or solidify the tetradecane, the COP of the chiller was 
worse and therefore the generation of cooling was more expensive. 
The tetradecane emulsion balanced out the additional generation 
costs with the reduction in distribution costs, resulting in the lowest 
specific cost for the three studied cases. 


7. Conclusions 


From this review it cannot be concluded whether PCM emul- 
sions and mPCM slurries improved the heat transfer phenomenon 
in comparison to water as a heat transfer fluid, as it depends on the 
combination of influential factors. The influence direction of these 
factors is known, but the degree of influence is still unknown. There 
is a lack of experimental studies, which are required to deepen the 
analysis carried out herein as well as to validate the results of dif- 
ferent numerical studies, as a great part of the studies has been 
verified through the same reference or with results that did not 
consider phase change. 

After analysis of the works compiled herein, Table 7 summarizes 
the objective magnitudes at the time of selecting a PCM emulsion 
ora mPCM slurry as heat transfer fluid or thermal storage material. 
The table includes the factors that influence the objective magni- 
tudes and how this influence occurs. 

This review has focused, on one hand, on the elaboration and 
fabrication of these new liquids, and on the compilation of the PCM 
slurries elaborated and utilized up to date in literature. On the other 
hand, this review focused on the analysis of properties, which is a 
key factor in the optimal development of storage fluids. The main 
problems in thermal storage material applications are subcooling 
and instability issues. In heat transfer fluid applications, it is inter- 
esting to achieve a higher heat transfer in comparison to a typical 
heat transfer fluid, such as water. Nevertheless, after systematic 
reviews, a clear conclusion still cannot be drawn, as several factors 
play a role. 

For the moment there are very few installation examples that 
make use of PCM emulsions and mPCM slurries. The main appli- 
cation seems to be oriented towards storage tanks combined with 
chilled ceilings, shifting the cooling production to the night in order 
to shave the demand peaks, improve COP, and make better use of 
nocturnal tariffs. There is also a study on the use of these new fluids 
in supply networks. 
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